Optical properties of quantum wires: 
Disorder— scattering in the Lloyd— model 
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The Lloyd model is extended to the exciton problem in quasi one-dimensional structures to 
study the interplay between the Coulomb attraction and disorder scattering. Within this model the 
averaging and resummation of the locator series can be performed analytically. As an application, 
the optical absorption in quantum box wires is investigated. Without electron-hole interaction, 
fluctuations in the well-width lead to an asymmetric broadening of the minibands with respect to 
the lower and upper band-edges. 



I. INTRODUCTION 

The optical properties of semiconductors near the 
band-edge are substantially influenced by the attractive 
electron-hole interaction which leads to pronounced ex- 
citonic lines and enhanced absorption above the band 
gap [Q . These phenomena become even more pronounced 
in systems of reduced effective dimensionality such as 
quantum-wells, wires or dots. In addition, the influence 
of disorder and impurity scattering increases with de- 
creasing dimension. The latter phenomenon, however, is 
considered usually phenomenologically by replacing the 
light frequency u; by w + 17, or convoluting the spectrum 
with an appropriate smoothing function of width 7. 

In this article, we investigate the interplay between 
excitonic effects and disorder-scattering in quasi one- 
dimensional structures, in particular quantum wires 
which possess a periodic step-like modulation along 
the one-dimensional direction, (quantum box structures, 
QBS), Fig. 1. Such structures have been proposed as 
favorite candidates to generate Bloch oscillations [§-Q, 
other examples are quasi one-dimensional molecular 
crystals [^-^|, yet no detailed theoretical study of the op- 
tical properties has been published. The optical proper- 
ties of one-dimensional systems are exceptional as even a 
small attractive electron-hole interaction leads to a dras- 
tic reduction in the oscillator strength at the gap energy 
so that the identification of the gap energy from optical 
absorption data is nontrivial. 

Our paper is organized as follows: Sect. II summarizes 
the basic description of the optical susceptibility in terms 
of the resolvent operator. In Sects. Ill and IV the av- 
eraging on the disordered chain is performed. Within a 
Wannier representation the averaging and resummation 
of the locator expansion can be done exactly, provided 
(a) only the lowest electron/hole subbands are consid- 
ered, and (b) a Lorentzian probability distribution for 
the disorder averaging is used. The latter is known as 
the Lloyd-model [10| which has been first used in con- 
nection with the calculation of the density of states in a 
disorded metal. Due to the simplicity in calculating var- 
ious averages this model has been found useful by many 



other researchers, in particular Refs. [[ri|~p5[. However, 
apart from Hoshino's work [Q on the dc conductivity 
(CPA, neglecting vertex corrections) all previous appli- 
cations of the Lloyd-model are exclusively restricted to 
single-particle quantities. Our work extends the Lloyd- 
model to a two-particle quantity and we give an exact 
solution for the optical absorption within a two-band 
model. Sect. V gives some analytical and numerical re- 
sults for QBS near the fundamental gap and Sect. VI 
contains our conclusions. 
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FIG. 1. Quasi-one-dimensional quantum-box-structure 
(left) and sketch of lowest valence/conduction subbands 
(right). (Not to scale). 



II. BASICS 

Our starting point is the excitonic contribution to the 
optical suszeptibility x{^) i^ ^ two-band effective mass 
approximation il|,|2f: 



xi^) 



2|dc- 



eoQ 



■E 



1 



huj + tS — En 



(fr'^nir,r) 



(1) 



where the nonresonant part has been neglected, (w > 0). 
ri is the volume of the sample, and |dcuP is the dipole ma- 
trix element between the valence and conduction band. 
^„(re,r;i). En respectively denote the electron-hole en- 
velope function and energy, and n labels discrete as well 
as the continuous part of the excitonic spectrum. 

[H^ + if, -I- V{r^ - rO] *n(n., rO = S„*„(r„, n.). (2) 



Hj, j=e,h denote the electron/hole Hamiltonians and 
V{r) is the attractive electron-hole Coulomb poten- 
tial. We assume a cylindrical wire of radius R with 
an infinite confining potential and restrict ourselves to 
the lowest electron/hole mini-bands with wave functions 
'^ jif) — u{r±)tpj{z). u{rj_) is radialsymmetric, nodeless 
and tpjiz) denote the wave function along the wire z-axis 
which obeys 



H, 



2mj dz^ 



v,{^). 



(3) 

(4) 



m- 



are the effective electron/hole masses and Vj{z) in- 
cludes the band offset as well as the localization energy 
due to the radial confinement. For a homogeneous wire 
the energy spectrum consists of free electron like mini 
bands with parabolic dispersions. These bands will be 
additionally split into subbands by the periodic modula- 
tion in composition along the wire axis. In the follow- 
ing we restrict ourselves to the lowest electron and hole 
subbands which are conveniently described in terms of 
Wannier-f unctions aj{z — pb) where p labels sites and b 
is the period of the box structure, Fig. 1. By an appropri- 
ate choice of the phases of the Bloch-functions 4'j{k,z) 
the Wannier-functions can be constructed in such a way 
that they are real, symmetric, and exponentially local- 
ized Ig^. Furthermore, the Wannier-functions of the 
lowest hole and electron subbands both have the same 
even parity. 

"^niZa, 2h)= ^ V'p"ph"<=(^<= ~ PM<^h(Zh - Phb)^ (5) 



Po,Ph 






(6) 



Within a next neighbour approximation the the matrix 
elements of the exciton Hamiltonian H = H" + TC" + V 
read 



Spa 



.Ph,Pc,Ph 



(E) 



n.k 



Ptp PIP p,p ' 

jjU) _ r (j) 
p,p ~ P>P I 

^P,P ^ *1 ■f'^P.P+l + ^P:P+l}- 



V. 



PcPhPcPh' (') 

(8) 

(9) 

(10) 



£0,4"' respectively denote the on-site energies and 
transfer elements between neighbouring wells and 
Vp^p^pji^^ « V{p, - Ph)Sp^pJpi,pt, is the electron-hole 
Coulomb matrix element. 

To calculate the optical absorption we rewrite Eq. (|l|) 
in terms of the resolvent operator G{E) 



x{tii^) « ~ yj Gp.,p.p,p (fi^ + iS) , 



p,p 



g{E) = (e-h 



(11) 



(12) 



where E ^ E + id, S > 0. 

Without disorder, the electron-hole wave function 
obeys the Bloch theorem and, hence, it can be labelled 
by a wave number k = 2ttk/N, k = 0, ±1, . . . ± N/2 






PiL + V 



^Po, 



Ph- 



(13) 



V — 1,2 ... N, where N is the number of sites on the 
chain. Periodic boundary conditions are implied. Hence, 



E — Enh 



Spc 



-i^,Ph+i^,Pj,+i/,Pi,+' 



^{E). 



(14) 



Furthermore, only fc = states V'p^'ph — '^p.-p contribute to the optical absorption 

.*(") i(n) 

J2gp,p,p,p{E) =nY,Gp,pME) =nY. ~ ° = NGooiE). (15) 

p,p p n "' 

The k — states are translational invariant and effectively describe a single particle in an external potential 

(ii + tl) [4% + 4-\) + (eo + eo) 4"^ + Vr^"^^ = E,,cl>i-^ , reZ. (16) 



Gr,f{E) is the corresponding Green-function which 
obeys the Dyson-equation 



G = Go + GnVG. 



(17) 



Boldface symbols denote matrices. Go is the Green- 
function for Vr = Q, and V is a diagonal matrix, Vrr> = 
VrSrr'- The formal solution of Eq. (17) is obtained by 



matrix inversion 



G = (l-GoV)-iGo 



(18) 



In addition to the attractive electron-hole interaction 
we shall consider fiuctuations in the on-site energies eg 
as well as in the transfer elements ti between adjacent 
wells of the Kronig-Penny potential. Such fluctuations 



arise from compositional and structural disorder of the radial symmetry so that mixing of higher subbands will 
wire and we assume that these fluctuations preserve the not be important. 

III. BARRIER FLUCTUATIONS 

Fluctuation in the barrier thickness lead to coupled fluctuations in the on-site energies and the transfer elements 
between sites q and q + I (nondiagonal disorder). In particular, we consider a correlated linear change of the onsite 
energies as mentioned by John and Schreiber [n2|. For the lowest subband ii is negative, thus, corrections to the onsite 
energies and hopping element have the same sign, i.e. Sq = aAt, a > 0. (The coupling of diagonal and non-diagonal 
disorder will be essential to guarantee the convergence of the locator expansion and non-negative absorption). For 
simplicity, fluctuations in the electron and hole parameters are assumed to be independent so that the Hamiltonians 
Eq. (0) are replaced by 

pS = 4'^ + ^^'^^i% + ^^'^^ii-i: (19) 

Ki = ^p^p+^ {4^' + Ki+i} + ^.P+i {4^'^ + A^s+i} ■ (20) 

In matrix notation H = T> + M +V, where P, A/" include electron and hole parts. To perform the disorder averaging 
the locator expansion of the resolvent will be used |2q] 

oo 

g{E) = Y,G^''HE), (21) 



^(»)(^) = _1 



n=0 

1 



E -V 



M-^ 



E~V 



^ ~ ■^p^""''pr""po""''p;""" ^ " '^p.p^pji^ 



Next, the disorder averaging will be performed term by term. Concerning the factors 



'PoPhPoPh 



PoPhPoPh 



four cases must be distinguished: 
1: (Po ^P,,Ph =pj 



■^PoPhPcPh _ '''^PoPhPoPh 



2: {Pa ^P,,Ph ^Pt,) 

■^poPhPoPh ''^PoPhPoPh + f^Ph.Ph+i {^1 + A^ph,ph+i } + '^Ph,Ph+i yi + A^Pi,,Ph+i I 



(22) 



where Ai = N + V. As I? is a diagonal matrix, we have {n > 2): 

^PcPhPeP. (^) - -f-f: ' (23) 

^ ■'^PoPhPoPh 

^W (E) = -. ^ ^^P^P^P.P. ^ (24) 

PcPhPcPh ^ ' Tp -T-) Tp •T-)_ ^ ^ 

-^ -^PoPhPcPh -^ -^PoPhPoPh 

, , , . . . 1 M (1) (1) M (1) (1) (2) (2) 

^("^ __(£;) = y^ ••• 'V PcPhPo 'pi Po 'pI 'p. 'pI 

11 I -<\ r -1-, E — 2?j, „,„ rn, E — 2? (1) (1) (1) (1) E — 2? (2) (2) (2) (2) 

pIpI Po"" 'pJ,"" ' PoPhPcPh p^^ 'pj^ 'p^ 'p'^ j,^ 'jjj^ 'pi 'pj^ 

M. (r.-2) (n-2) (ii-l) (n-1) Ai (n-1) (n-1) 

Pe Ph Pe Ph pi Ph PePh 



(25) 



(26) 



(27) 



(28) 



3: {Pa ^P,,Pi. ==Ph) 



M 



PoPhPaPh 



4: [p^ ^p^^ pi^ ^ pjj analogous to case 1. 

The averaging procedure requires the evaluation of the following multiple integrals 



^PePhPcPh ^ 






dAtl^ . . . dAt"^i dAtj2 • • ■ rfAi^i P(Atl2, r) • ■ ■ P{At%^,r) X 
X P(Ai'i2, 7") . . . P{At-^,,j-)gl'^l--^^{E, Atl„ . . . Afk,,At\^, . . . Ai^^), 



(29) 



(30) 



where P(A, 7) is the probability distribution function 
on the disorder configurations. In the Lloyd-model this 
function is a Lorentzian 



7 



^(A,7) = - 



1 



TT A2 + 72 



(31) 



where 7 parametrizes the width of the distribution. 

The integrations in Eq. (^0|) will be performed step by 
step in the indicated sequence. In Eqs. (pS-Ea) each fac- 
tor is a holomorphic function of Atifj in the lower complex 
plane including the real axis, in addition, this function is 
bounded if a'" ^ 0. Thus, the averaging can be performed 
by contour integration closing the integration path in the 
lower Ai'j'g-plane where Eq. IBWj has a first order pole. 
This just leads to a replacement of At!j'2 by —ij'". Fur- 
ther integrations follow the same reasoning. Hence, the 
resummation of the locator series is trivially possible and, 
as a result, we obtain 



g:={g{E)) = (E-H 



(32) 



TC is an effective Hamiltonian where the electron/hole 
parts of Eqs. (0) are replaced by 



7^S = '5p.p(e?^-2z7(^^«) 



(^P,P+i + %p+i) (4'^ - *7^^'^) 



(33) 



Thus, the influence of barrier-disorder is captured by 
the resolvent of the pure wire where e^'jil ■"■" ^^ 
placed by complex parameters e 



U) 



Aj) 



t 



U) 



Aj) _ 2i7(j)a(j'), 
i^U) ^ Although TC is no longer hermitian 
eigenstates exist and, moreover, these are of the same 
form as for the pure wire, in particular the eigenvalues 
are 



£;0")(fc) = ei^') + 2t[^^ cos{kb). 



(34) 



For t^ < the lower/upper band edges are at kb = O/tt 
so that the complex electron/hole energies are distributed 
between e„-2|ii|- 217(0 + 1) and e^ + 2\ti\-2t'y{a-l). 
The parameter a must be restricted to a > 1 for both 
electrons and holes, otherwise there will be energies with 
a positive imaginary part, i.e. the Green-function will 
have poles in the upper energy plane which would violate 
causality. In addition, the optical absorption would be- 
come negative in some part of the spectrum. Note, that 
the correlated fluctuations in the diagonal/nondiagonal 
disorder will not simply lead to a Lorentzian smoothen- 
ing of the absorption spectrum. In particular, for a — 1 
the correlated shift of the on-site energies and transfer el- 
ements exactly cancel at the upper band edge so that the 
square root singularity of the density of states survives if 
electron-hole interaction is omitted. 



IV. WELL FLUCTUATIONS 



Fluctuations in the well sections predominantly lead to fluctuations in the onsite energies e^. For convenience, 
fluctuations in the hopping elements will be omitted, however, a correlation of the electron/hole energies will be 
retained: A^ = aA", a « m^/mh > (diagonal disorder). 

To perform the disorder averaging a regrouping of the Hamitonian as 7Y = S + A^ is necessary, where 



■'^PePhPoPh ~ "PhPh ("PoPo^o + ^1 ("Po,Pc + l + "Pe,Pc-lj| 

"I" ^PbP" ("PhPh^o + ''I (,"Ph>Ph + l + "Ph.Ph-lj J '^ '^PePhPePh' 
"PoPhPoPh ~ "PhPhApjj + "PePoAp^. 



(35) 
(36) 



Using the locator expansion Eqs. (ElU22) with V replaced by B the matrix elements of the resolvent become (n > 2) 



c(°^ s - 

^PoPhPoPh 



Po'Pi.Po.Ph 



■^PoPhPePh 



1 



M 



PcPhPePh 



£; - (Ap^ + aApJ S - (Ap^ + aApJ 

1 



(37) 
(38) 



E-- E 



•^PoPhPcPi 



■^pJpiPoPh 



_^ E - {Ap^ + aApJ E - {Ap 



- aA„i ) E 

Ph' 



(Ap 



tA 



Ph^ 



A4 n-2 n-2 „-l „ 
Pc Ph Pc P., 



p? Vh VcPh 



(«) _ _ 

PcPhPcPh 



(^,Ai 



,An,). 



(39) 



(40) 



^ - (Apr^ + "^prO E - (Ap^ + aApJ 
The disorder averaging of Q is done along the same route as before 

/CO 
dAi •••dAjvP(Ai,7) •••P(A^,7)a 
-oo 

As a result, we obtain 

g{E) := (giE)^ ^{E + i{l + a)7 - M)~^ . (41) 

Thus, with respect to the pure system, diagonal disorder just leads to the replacement oi E = huj by the complex 
quantity E ~^ E + i{l + a)7, which leads to the expected Lorentzian broadening of the optical absorption spectrum. 



V. APPLICATIONS 

The Green-function of the noninteracting pure wire 
is well known from textbooks, e.g. Economou |2^]. In 
a sligthly rewritten form which is suitable for the ana- 
lytic continuation to complex Tq — el + e|J, Ti = i\ + i'l, 
i^Ti < 0), it reads 



cm) 



pi'-"i(^) 



[E - To)2 - 4T2 



(42) 



where uj' ^ [hu ~ To)/(2|Ti|), vq = \Vo/{2Ti)\. W con- 
taines the interband dipole matrix element and an over- 
lap integral between Wannier- functions. A similar result 
has been obtained by Ishida et al. [0. Note, kuj is mea- 
sured with respect to the gap energy Eg, of the homoge- 
neous wire. Fig. 1. 

For a quantum wire we have to replace the three- 
dimensional Coulomb potential by the envelope averaged 
potential which can be approximated quite well by [y 



where 



V{z) = - 



1 



^iTeoe\z\+(3R 



(47) 



Pi{E) ={E-To-^{E- To)2 - ATlj /(2Ti). (43) 

\/x denotes the square root whose imaginary part has 
the same sign as Qx, see Eqs. ( ^ , pl| ). 

First, we consider a local electron-hole interaction, 
Vr — VoSr,o, Vq < 0. In this case Eq. ( p7[ ) can be solved 
analytically 



GooiE) 



GUE) 



l-VoG"„o{E) 



(44) 



Without disorder (i.e. real To,Ti) the optical absorption 
X2 = '^X of ^ wire becomes 



X2{^) 



\w\ 



^ 






W\ < 1, 



2|r,i yr 



(45) 
(46) 



where R is the wire radius, /3 « 0.3, and e is an aver- 
age dielectric constant. For a GaAs based wire of radius 
R — 50A, V{0) « — O.leV which is approximately the full 
electron-hole transition bandwidth 4|Ti|, hence vq « 0.5. 
To solve the Dyson equation Eq. (n^) the effective ex- 
citon potential Eq. ( p7| ) is truncated to a finite range, 
Vr — V{rb) — 0, |r| > s. Then, the optical absorption, 
Eq.(n2|), is obtained numerically by solving a (2s + 1) di- 
mensional linear set of equations for Gio{E), I = —s, . . .s. 
Figs. 2-6 display various examples for the perfect and dis- 
ordered wire. 



VI. CONCLUSIONS 

The optical absorption in a quasi onedimensional sys- 
tem sensitivly depends on both the electron-hole inter- 
action and the disorder scattering. Diagonal disorder 
just results in a Lorentzian smearing of the optical lines. 



whereas correlated diagonal/ofFdiagonal disorder is differ- 
ent: the broadening is stronger at the lower than at the 
upper band edge, Fig. 2. However, the shape of the ab- 
sorption spectrum is drastically changed even by a tiny 
electron-hole interaction: The divergences at the band 
edges are completely removed so that the identification of 
the gap is a nontrivial problem in quasi one-dimensional 
semiconductors. Figs. 3-5. For a local electron-hole in- 
teraction there is only but a single exciton bound state. 
Fig. 3, whereas the (truncated) Coulomb potential leads 
to an additional pronounced spike near the gap and an 
oscillatory structure within the absorption band. Fig. 6. 
These structures are insensitive with respect to the type 
of disorder. 




FIG. 4. Optical absorption with varying electron-hole 
interaction. Fixed barrier-disorder, q° = a'' — 1.1, 
7° = 7*' = 0.025. Energy scales as in Fig. 2. 




-2-1 1 

CO' 

FIG. 2. Optical absorption of the noninteracting elec- 
tron-hole system with varying barrier-disorder. 
a" — a^ = 1.1. Energies are given in units of one half of the 
transtion band width, lower band edge is at lu' — —1. 




FIG. 5. Optical absorption with varying electron-hole in- 
teraction. Fixed well-disorder, (1 -|- 0)7 — 0.1. Energy scales 
as in Fig. 2. 
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FIG. 3. Optical absorption of the pure wire with variation 
of the local electron-hole interaction. Energy scales as in 
Fig. 2. 
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FIG. 6. Optical absorption with varying range s of the 
truncated electron-hole Coulomb interaction, wo = 1. Fixed 
well-disorder. Q° — a^ = 1.25, 7° = 7'' = 0.005. Energy scales 
as in Fig. 2. 



Our approach is based on the locator expansion of the 
Green-function (v. Neumann series) whose convergence 
is rarely established in general. A more rigorous treat- 
ment, however, shows that Eq.(B^ is indeed correct if 
a > 1 1^. For < a < 1, Eq.p3|) leads to a negative 
absorption in some part of the spectrum which is defini- 
tively incorrect. It seems that the Lloyd model cannot be 
formulated for purely nondiagonal disorder, i.e. a = 0. 
A similar problem is known for the extension of the CPA 
to nondiagonal disorder [bI| , but it may be also due 
to the pathological Lorentzian probability distribution, 
Eq. (|l|), of the Lloyd-model. 
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